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I. INTRODUCTION 

It is well known at present that the polaron problem has much more broad significance than simply a model for 
the partial form of interaction between electron and phonons in the ionic crystal as it was introduced by Frohlich pj . 
Hamiltonian of the polaron is considered mostly as a fundamental model for interaction between particle and quantum 
field where various non-perturbative methods for the quantum field theory can be verified in the entire range of the 
coupling constant a. 

Like any other quantum system the polaron can be described both on the basis of the solution of the Schrodinger 
equation and by means of calculation of the corresponding Feynman functional integral. The first way allowed to 
introduce the idea of the self-localized polaron and to find the exact asymptotic value for the ground state energy 
Eq of the system in the limit a — ► oo But the great advantage of the Feynman variational principle for the path 
integrals is the possibility to calculate Eo(a) as the continuous function for any a and to find the lowest estimation 
for the polaron binding energy in the intermediate coupling regime 4J. 

There were a lot of attempts to calculate the ground state energy by means of variational principle for the 
Schrodinger equation (usually it is called the operator calculation) and some of the trial functions led to pecu- 
liarity for the function E (a) near the point a ~ 7. These results caused the discussion about existance of the "phase 
transition" between two qualitatively different states of the polaron (see review Q on this problem). In the series 
of papers cited in [{| it was proved that the function Eq{ol) was analytical for any a and "phase transition" did not 
exist. But it is important to stress that this strict proof did not consider the constructive computational algorithm 
for analysis of the self-localized and another states of the polaron. 

In the present paper we found new variational ansatz for the ground state vector for the considered system which 
leads to the function Eo(a) fairly well coincided with Feynman's calculation. It seems to us that the results make the 
question about the ground states of the polaron much more clear and descriptive. 

It will be shown below that there are two collective states of the polaron system which describe the "dressed" 
electron and do not connect with excitation of the free phonons. In general case each of them is the mixture of 
delocalized and self-localized states of the "bare" electron. The energies of these states are splitted so that the both 
eigenvalues prove to be the continious functions of a and the ground state energy is similar to the result found by 
Feynman. In the limit a — > oo the corresponding wave function transforms eventually to the self-localized state. But 
it is important that both states exist in the entire range of a and a pseudo-intersection between these energy levels 
could be the reason for imaginary "phase transition" for some sort of variational estimations. 

Our analysis also showed that the conception of self-localized state proved to be applicable for more general models of 
the particle-field interaction even in the case when usual perturbation theory includes both the infrared and ultraviolet 
divergencies. 

II. DESCRIPTION OF TWO BASIC STATES 

We use the standard form of the Frohlich Hamiltonian (fi = m = uj = 1) 

H = ~\ A + £ aW + ^ (f ) V2 E I (<^X + ^ k X + ) ■ (D 

k k 

Here r is the coordinate of electron; a k and are the operators of annihilation and creation of the phonons with 
the momentum k; a is the dimensionless coupling constant. 

The problem is to find the common solution of the Schrodinger equation 

H\*)=E(a,P)\V), (2) 
and the equation corresponding to the conservation law of the total momentum 

P|*> = (-iV + 5>4a k )l*> = P|*>- (3) 

k 

Let us consider now two well known approximate variational solutions of these equations which we will use as a 
basis for our ansatz. We remind shortly the algorithms for obtaining mentioned solutions. 
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The first one is defined as follows |6J 

|^ ) (r)) = -^exp^Pr + ^z; k (P)( e -* r a +- e ^ ak )^ |0); a k |0) = 0. (4) 

We call this state delocalized (D) because its structure corresponds to the polaron where the cloud of the virtual 
phonons follows the freely moving electron. The state is the eigenfunction for the operator of total momentum and 
the parameters «k(P) are defined by the minimum condition of the energy found as the average of the Hamiltonian 
JTJ with respect to this state 

%(P) = -^(f) 1/2 t[fc2 _ 2(p '_ Q)k + 2r Q = X>J(P>. (5) 

For simplicity we'll consider all further calculations for the ground state with P = when the energy is defined by 
well known simple formula 

E D = H DD = J dv{^ D \v)\H\¥ Q D \v)) = -a. (6) 

The second basic state is generated by the product ansatz which was firstly introduced by Pekar Q, who realized 
Landau's [j| suggestion that an electron in an ionic crystal can be trapped by the static lattice distorsion produced 
by the electron itself. This product ansatz was proved to be the asymptotically exact solution of the Schrodinger 
equation in the limit a — ► oo as it was shown by Bogoliubov Q . 

We can write it as follows 



|*(r,R)) -0(r-R)exp^ Uk ( e - kR a+-e ikR a k )^ |0>; J dr|0(r)| 2 = 1. 



(7) 



Here R is the arbitrary point in the space where the electron is localized. The classical component of the phonon 
field Mk is also defined by the variation of the energy functional and has the following form ([3j) 

uu = -2^(^) 1/21 -Jdr\^e'^, (8) 

and the wave function of the electron in the ground state should be found as the solution of the following nonlinear 
equation with eigenvalue Ep 

h U-*l> a UttiH!^-^ /*, [ dr M*^'f -E pW) =0. (9) 
2 J |ri — r| J J |ri — r 2 | 

In spite of the exact solution of this equation was found numerically we will use for further calculation the 
Pekar's function 0o(r) which is a very good approximation 0,0 



- lT - , ^/ l4 + 42 8 C+ 45c 2 ( 1 + 6? - + CbV ^ ^ 



with 

b = 0.931307a, c = 0.451668. 
The optimal value for the variational parameter 6, c leads to the polaron binding energy 

E P = -0.108504a 2 
which is very closed to the accurate value 0,0 —0.108513a 2 . 
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It is clear that the set of the state vectors (JTJl correspond to the localized states (LS) of the electron which can 
be considered as the states with local breakdown of the translational symmetry of the system in the sense that was 
discussed in the fundamental Bogoliubov's paper {u%. The energy of the system should not depend on the point 
of the electron localization R because of this symmetry, and in order to reconstruct it let us use 

nulla 

the linear 

combination of degenerate states |W(r, R)) 



with 



|*W(r)) = ^= J dRe iPR |*(r,R)), 
Nl= j dR j dr(j){r) 4>{r + R) cxp[^ ^ ( e - 4kR - 1)]. 



(11) 



This state vector leads to the corresponding estimation for the polaron binding energy (we fulfilled real caculations 
for the case P = and used JTDJl) 

El =H LL = J dr{¥ L \r)\H\¥ L \r)) (12) 

Fig.l compares the functions Ed(o), El (a) with the value Ep(a) calculated with the state vector (JJJ without 
taking into account the translational symmetry. 

III. NEW ANSATZ FOR THE POLARON STATE 

Fig.l shows that generally speaking there are no obvious restrictions on the coupling constant a for both continuous 
functions Ed (a) and El (a) corresponding to two qualitatively different states. The only suspicious point nearby a ~ 7 
is defined by the intersection of these functions where the states are degenerated. The similar situation is known in 
the zone theory of electronic states in crystals as the intersection of resonant zones 0|. We use the same method in 
order to analyze the transformation between these states by means of the following linear combination 

|$ P (r))=C 1 |4 D) (r) + C 2 |4 i) (r)) (13) 

with only linear variational parameters and fixed values of internal parameters for the basic states. 

In a series of previous variational calculations of the polaron binding energy for the intermediate coupling constant 
various rather complicated trial functions for the localized or delocalized state were used separately (see, for example, 
the review ^^) but the simple state vector (|13Jl can't be reduced to any of them. 

The subsequent calculations are rather bulky but standard. By the definition the both basic vectors are eigenstates 
for the operator of total momentum and the ansatz i|13|) should be used only for the functional corresponding to the 
Schrodinger equation The condition of vanishing the determinant obtained after variation of this functional on 
the coefficients C\ 2 leads to the following formula for two eigenvalues: 

E ± {a) = 2(1 ^ g2) (Hdd + H LL - 2SH LD ± E R ) 



E R = V(Hdd + H LL - 2SH LD y - 4(1 - S 2 )(H DD H LL - H LD H LD )- (14) 

Here Hdl, Hdd and Hll are the matrix elements of the Hamiltonian between the corresponding basic states 
and S is the ovelapping integral between this states being nonorthogonal in general case. Omitting the calculations 
in detail we write out the final analytical results for these values when P = 0. 

The most simple value is 

H DD = -a. (15) 



The energy of localized state is defined by more complicated formula 
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FIG. 1: Energies of polaron for Pekar's (Ep), delocalized (Ed) and localized (El) states. 



with 



At last 



H LL = j^ JdRj dr«Mr + R)exp(^<4e- ikR ) 
(4 A + E ^ e " <kR + ^ (f ) V2 E ^ ^ + ^ tk(r+R) )}^o(r) 

k k 

N 2 = JdnJ dr</> o (r)0 o (r + R)exp(^u k e^ kR ). 



-a/4 



[ drcf> (r){-a + ^ (^ k?i k v k ) 2 } cxp(J^ u k i; k e jkr ) 



(16) 
(17) 



(18) 
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and 

5=—— / dr0 o (r)exp(^ Uk z; ke jkr ). (19) 

k 

Formulae (|14fl - (|19|l permitted us to calculate the functions E±(a) in the entire range of the coupling constant 
and Figure 2 shows the results. The well known estimation for the polaron energy Ef{o) found by Feynman jj] is 
drawn for comparison. One can see that the functions E- and Ef are very similar continuous functions. Certainly, 
_E_ becomes lower than Ep for rather large a because of more accurate main asymptotical coefficient. The value 
doesn't depend essentially on the detailed form of the electron function </>o(r) as one can see from the curve Eg 
corresponding to the Gaussian approximation for this function. 




FIG. 2: Polaron binding energies calculated by Feynman (Ef) and on the basis of our ansatz (E±); Eg is the analog of E- 
but with the Gaussian electron wave function. 

The noticeable difference between our result and Feynman approximation in the intermediate range is defined by 
more accurate constant term in the strong-coupling series and it can be improved by means of the second order 
approximation |15| . Figure 3 confirms it on the example of one-dimensional polaron with Hamiltonian 

^=-^ + E a K +2 1/4 (f) 1/2 £( e "X +e-^at) (20) 

k k 
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where the role of the constant term is essentially less and we used the exact solution of the nonlinear Schrodinger 
equation for localized state. 



E \ 

1 2 3 4 5 




FIG. 3: The same as on Fig. 2 but for the one-dimensional polaron. 



IV. THE MODEL WITH DIVERGENCIES IN PERTURBATION THEORY 



The polaron problem is the unique model because all diagrams of the quantum perturbation theory for it are finite 
without renormalisation procedure . Therefore it is of interest to apply our conception of localized states to more 
general system. As an example, we consider the following Hamiltonian 

= -^A + ^ Wfcfl +« k + -| E (e' kr « k + e- kr «+) , Wfc = fc=|k|. (21) 

k k 

It is easy to check that the particle self-energy diagram of the second order is infinite in this model because of 
ultraviolet divergence of the integral over k. As a consequence, the delocalized state analogous to @ does not exist 
here. 
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Let us now consider the localized states analogous to (lllfl 

10) (22) 



|*(r,R)) = 0(r-R)exp ^( 



-ikR + _ i 2 
U k e °k 2 k 



with 

and the wave function </>(r) satisfies the nonlinear equation of the same form @ as for polaron. 
The total energy of this states 

E = [dr #r)(-i A)#r) - ]T w fc i£ (23) 

k 

doesn't contain the ultraviolet divergence due to the particle formfactor and proves to be the finite value for any 
coupling constant g. 

But there is another problem here connected with infrared divergence. The matter is that the state vectors (1221 are 
not normalized because of the term ^ k u\ which is proportional to the total number of phonons and diverges over 
k — ► 0. The situation changes radically when we reconstruct the translational symmetry of the system by analogy 

with mi 



| * ^ (r)) = _L= JdR^v-K) exp ({PR + £ (■ 



"ke-^-kl I l°>- ( 24 ) 



2 

The integral of normalization in this case (for simplicity P = 0) is following 

N, 



I = JdRjdv 0(r)0(r + R) exp ^ u 2 u ( e -* R - 1) j (25) 



and includes only the converging integrals. 

Figure 4 shows the binding energy of the "dressed" particle 



E L {g) = J dr{^\v)\H\^\v)\, 

which is well defined and finite in the whole range of the coupling constant. At the same time this value can't 
be calculated on the basis of usual perturbation theory because of the same reason as for the bound state in any 
potential. 



V. CONCLUSIONS 



The main purpose of this article is to make clear the question of "phase transition" which appears for a number 
variational approach used for description of the polaron in the Schrodinger picture of the problem. Actually this 
question includes two sides 

(1) Do qualitatively different polaronic states really exist? (2) Does the jump-like transition between these states 
take place for some value of the electron-phonon coupling constant? 

The new variational ansatz for the Schrodinger was found in our paper and it was shown that the energy of the 
polaron ground state Eq (a) is continuous function of the coupling constant in accordance with the analysis on the 
basis of the Feynman path integrals. But at the same time there are two qualitatively different states of the polaron 
without excitation of real phonons. The pseudo-intersection between these states describes in detail the transition 
between the delocalized state of the electron in crystal for small a to the localized one in the strong coupling limit 
and can lead to some physical effects nearby the point of quasi-degeneracy. 

One of the advantages of the considered operator calculation of the value Eq{q) is defined by the fact that we 
use only the eigenfunctions of the total momentum operator and it permits us to calculate the effective mass of the 
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FIG. 4: Self-energy El of the "dressed" particle for the model with the Hamiltonian (21) (Ep is the energy without taking 
into account the translational symmetry). 

polaron by consistent way. We'll consider the corresponding rather bulky calculation separately Besides, the operator 
approach allows to find the exact asymptotic serieses in the both limit cases on the coupling constant and calculate 
by regular way the subsequent corrections in the intermediate range. 

And the important result of the paper is defined by the generalization of the conception of localized states in 
application to the system with the ultraviolet and infrared divergencies of the perturbation theory. We beleive that 
the considered non-perturbative description of the "dressed" particle is of great interest for quantum field theory. 
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